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Abstract. Big Bang Nucleosynthesis (BBN), along with the cosmic 
background radiation and the Hubble expansion, is one of the pillars 
of the standard, hot, big bang cosmology since the primordial synthe- 
sis of the light nuclides (D, ^He, ^He, ''Li) must have occurred during 
the early evolution of a universe described by this model. The overall 
consistency between the predicted and observed abundances of the light 
nuclides, each of which spans a range of some nine orders of magnitude, 
provides impressive support for the standard models of cosmology and 
particle physics. Here, the results of recent, statistically consistent tests 
of BBN are described. This new confrontation between theory and data 
challenges the standard model. The crises confronting BBN are identified 
and several possible resolutions are outlined. 



1. Introduction 

The discovery of the Cosmic Background Radiation (CBR) by Penzias &: Wilson 
(1965) transformed forever the study of Cosmology from an exercise in philoso- 
phy to the pursuit of science. The presence of the CBR in an expanding Universe 
favors the hot big bang cosmology. A Universe described by this model was very 
hot and very dense during early epochs in its evolution. As a consequence, it is a 
prediction of this "standard" cosmological model that, briefly, the early Universe 
was a primordial nuclear reactor in which the light nuclides D, ^He, ^He and ^Li 
were synthesized in astrophysically interesting abundances (for details and ref- 
erences see, e.g., Boesgaard & Steigman 1985; Walker et al. 1991). Thus, along 
with the CBR and the "Hubble" expansion. Big Bang Nucleosynthesis (BBN) 
provides one of the three pillars supporting the standard model of cosmology. 
The standard hot big bang model is, in principle, falsifiable. In contrast to cos- 
mology as theology, this empirical model is not a matter of faith but, rather, 
demands our eternal vigilance and critical scrutiny. The success of the model is 
gauged by the degree to which the BBN predictions are consistent with the pri- 
mordial abundances of the light nuclides inferred from observational data. Over 
the years BBN has emerged unscathed from the confrontation between theory 
and observations, providing strong support for the standard, hot big bang cos- 
mological model (e.g., Yang et al. 1984; Boesgaard & Steigman 1985; Walker et 
al. 1991). This success has, however, not spawned complacency, and the testing 
continues. In recent years, as the astronomical data has become more precise, 
hints of a possible crisis have emerged (Copi, Schramm & Turner 1995; Olive &: 
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Steigman 1995; Hata et al. 1995). It is my goal here to describe the impressive 
success of BBN and to map out the paths leading to the current challenges to 
the standard model. To better appreciate these challenges and the opportunities 
they present, an historical analogy may be instructive. 

1.1. Three Crises For The 19th Century Standard Model 

The gravitational theory described by Newton (1686) was outstandingly success- 
ful in explaining the motion of the moon and planets. Perhaps one of the most 
thoroughly tested physical theories in history, Newtonian gravity had become 
the standard model of 19th Century Physics. Soon, however, some challenges to 
the standard model emerged. The nature of these challenges and their different 
resolutions provide some interesting lessons for the emerging crisis in BBN. 

(i) Perturbations to the Orbit of Uranus 

Deviations in the orbit of Uranus from the predictions of Newtonian gravity 
(the standard model) led Adams and LeVerrier to predict the existence and 
location of Neptune. The standard model was used to discover something new, 
verifying the accuracy of the data (the orbit of Uranus) and providing spectacular 
support for Newtonian gravity. 

(ii) Perturbations to the Orbit of Neptune 

Observations of the newly discovered Neptune suggested that its orbit, too, 
was being perturbed away from the standard model predictions. So began the 
long search which culminated in the discovery of Pluto. Pluto, however, is 
not responsible for measurable perturbations to the orbit of Neptune - it is too 
small. Rather, here we have a case of insufficiently accurate data. The discovery 
of Pluto was serendipitous; more accurate observations of Neptune's orbit are 
entirely consistent with the predictions of Newtonian gravity. 

(iii) Precession of the Perihelia of Mercury 

By the mid-19th century LeVerrier had noted a discrepancy between the 
predicted and observed precession of the perihelia of Mercury. LeVerrier and 
others proposed one or more planets (Vulcan) between Mercury and the Sun to 
resolve this crisis. None were found. Alternately, it was proposed by Newcomb 
and others that the perturbing mass might be in a ring of dust or asteroids. 
This, however, would have perturbed the orbits of Mercury and Venus in conflict 
with observational data. A more radical solution, modifying the inverse square 
law, was proposed by Newcomb (1895). This, however, is in conflict with the 
accurately observed lunar orbit. 

As is so well known, the resolution of this crisis confronting the 19th century 
standard model was new physics! Einstein's General Theory of Relativity (1916) 
predicts a precession in beautiful agreement with that observed. 

Three crises, three different resolutions: the standard model preserved and anew 
discovery (Neptune); the standard model preserved and insufiicently accurate 
data (Neptune/Pluto); the standard model replaced (perihelia of Mercury). 
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2. Consistency of the Standard Model 



2.1. Predictions 

Now let us turn to the standard model of cosmology and the predictions of 

primordial nucleosynthesis. Employing measured weak interaction rates and 
nuclear reaction cross sections the primordial abundances of the light nuclides 
are predicted by BBN as a function of only one adjustable parameter, r/, the 
universal ratio of nuclcons (baryons) to photons [t] = Nb/Nj;7]iq = lO^^r)). 
The predicted abundances of ^He {Y is the ^He mass fraction) D and ^Li (2/2 = 
Nd/NhjVj = Nu/Nh) are shown for 1 < ??io < 10 in Figure 1 from Hata et al. 
(1995). For clarity of presentation the predicted abundance of ^He, very similar 
to that of D, is not shown. 

The predicted abundances depend on the universal expansion rate, t^^, 
during the epoch of BBN (~ 3MeV> Tbbn^ 30A;eF;0.1< tsBN^ lO^sec). For 

the early Universe oc p^^Qj-, where pxoT is the total mass-energy density. For 
the "standard" model (SBBN), proT is dominated by photons, electron-positron 
pairs and three flavors of light, left-handed neutrinos (fg, f^, J^r). 



In (1), p^ is the contribution from one flavor of light [m^ << Tbbn) neutri- 
nos. To account for a possibly massive r-neutrino and/or for other, new particles 
beyond the standard model, it is convenient to modify eq. (1) by introducing 
N^, the "effective" number of equivalent light neutrinos (Steigman, Schramm &; 



For SBBN, Ny = 3; for Ni^ 7^ 3 the universal expansion rate at BBN is modi- 
fied. For Ny > 3, the universe expands more rapidly leaving less time for the 
conversion of neutrons to protons. Since most neutrons are incorporated in ^He, 
Ybbn increases with N^, (and, vice- versa). Therefore, it is convenient to use Ni, 
as a second parameter to explore deviations from SBBN and extensions of the 
standard model of particle physics (Steigman, Schramm &; Gunn 1977). The 
results in Figure 1 are for SBBN {N,y = 3). 

For 1 < r/10 < 10, the predicted abundances of the light nuclides span 
a range of some 9 orders of magnitude from ~ 10"^*^ — 10"^ for Li/H, to ~ 
10-5 - 10-^ for D/H and =^He/H, to ~ 0.1 for ^He/H. 

2.2. Observations 

Primordial abundances are, of course, not observed. Rather, they are inferred 

from astronomical data. Some, such as D and ^He, have been mainly observed 
"here and now" (in the solar system and the interstellar medium (ISM) of our 
own Galaxy). For these nuclides it is necessary to extrapolate from here and 
now to "there and then" to derive their universal primordial abundances. ^Hc 
and ^Li are observed (in addition to here and now) in regions where much less 
chemical processing has occurred (low metallicity, extragalactic HII regions for 
^He; very metal-poor halo stars for ^Li). For these nuclides the extrapolations 
to primordial abundances are smaller. 




(1) 



Gunn 1977). 



PtOT =Pj + Pe + N,p\ 
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Figure 1. The SBBN predicted abundances (solid lines) of ^He (Y 
is the ■*He mass fraction), D (7/2 = D/H), and ''Li (yj = Li/H) as a 
function of the nucleon-to-photon ratio rj. The dashed lines are the 
la theoretical uncertainties from Monte Carlos. The shaded (dashed) 
contours are the regions constrained by the observation at the 68% 
(95%) CL. 



In addition to observational uncertainties and those associated with the 
extrapolations to primordial abundances, systematic effects in deriving abun- 
dances from data may contribute to the overall uncertainties. The bad news 
is that such systematic uncertainties are difficult to constrain. The good news 
is that the sources of possible systematic errors are different for the different 
nuclides. 

2.3. Testing SBBN 

The relatively strong and monotonic 1/2 vs. rj relation visible in Figure 1 points 
to D (and, to a lesser extent, ^He) as an ideal baryometer. If the primordial 
abundance of D were known, for example, to ~ 40%, the universal density of 
baryons would be known to ~ 25%. The large extrapolation from here and now 
to there and then has inhibited the implementation of this approach. Rather, to 
avoid this large extrapolation, a more conservative approach has been adopted. 
Since any D incorporated in a star is burned (to ^He) and there are no significant 
astrophysical sources of post-BBN D, the abundance of D observed anywhere at 
anytime provides a lower bound to its primordial abundance (e.g., y2P > 2/20; 
y2P > y2iSM)- From Figure 1 it is clear that a lower bound to y2P leads to an 
upper bound to rj. 

It is difficult to avoid the uncertainties of chemical evolution models in using 
observations of D to infer an upper bound to y2P- However, Yang et al. (1984) 
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noted that since D is burned to ■'He and some ^He survives stellar processing, 
the primordial abundances of D + ^He are strongly correlated with the evolved 
abundances of D + ^He. Burying the stellar and evolution model uncertainties 
in one parameter, 53, the "'He survival fraction, Yang et al. (1984; also. Walker 
et al. 1991) used solar system data to place an upper bound on primordial D 
(and/or on D + ^He). An upper bound on y2p provides a lower bound on 77 (see 
Fig. 1). 

Due to the "valley" shape in the BBN prediction of Li vs. r] (sec Fig. 1), an 
upper bound to y-jp will provide both lower and upper bounds to r]. The lithium 
abundance also offers a key test of the standard model since its primordial value 
must not lie below the minimum predicted {yjBBN^ 1 x lO"^'^). 

One test of the consistency of SBBN is to use D, ''He and '^Li to infer lower 
and upper bounds to rj {tjmin, 'Hmax) and to check that rjMiN < Vmax- If 
SBBN passes this test, the "^Hc test" may be applied. The predicted ^He mass 
fraction, Ybbn, is a very weak function of rj, increasing from Ybbn = 0.22 at 
7710 = 1 to Ybbtv = 0.25 at 7710 = 10. Thus, it is key to the success of SBBN 
[Ni, = 3) that for tjmin < f] < rjMAX, Yp (the inferred primordial abundance) 
is consistent with Ybbn{ii) (the predicted abundance). 

2.4. Consistency 

Yang et al. (1984) were among the first to carry out a detailed analysis of the ob- 
servational data and to implement the tests described above. From D, ^He and 
^Li (with gs > 1/4, see Dearborn, Schramm &; Steigman 1986) they found con- 
sistency: 3^ ??io~ 7, leading to a predicted range for ^He: 0.24^ Ybbn^ 0.26. 
Comparing with the rather sparse data available, they derived 0.23^ Yp^ 0.25 
and concluded that SBBN passed the ^He test. They did note that SBBN is, in 
principle, falsifiable and pointed out that if future comparisons should increase 
Vmin and/or decrease Yp, consistency would require < 3, modifying the 
standard model. 

By 1991 uncertainties in the neutron lifetime (as well as its central value) 
had been reduced considerably permitting a very accurate prediction of Ybbn 
vs. r] (at the 2a level, Ybbn is known to it 0.001; see Thomas et al. 1995). 
At the same time there was extensive new data on lithium (in halo stars) and 
helium-4 (in extragalactic HIT regions). Applying the above tests. Walker et al. 
(1991) found 2.8 < 7710 < 4.0 and 0.236 < Ybbn < 0.243. From the HII region 
data. Walker et al. (1991) derived Yp = 0.23 ± 0.01 and concluded that SBBN 
passed the ^He test. However, they did emphasize, "that if our lower bound on 
T] were increased from rjio = 2.8 to 7710 = 4.0, the window on Ni, would be closed 
(for Yp< 0.240)." 

2.5. Crisis? 

Recent applications of the two consistency tests {tjmin < Vmax'^ Yp = Ybbn'^) 
have provided hints of a possible crisis (Copi, Schramm & Turner 1995; Olive & 
Steigman 1995; Hata et al. 1995). The two "weak links" are the lower bound 
on T] inferred from D and ^He observations and the upper bound on Yp derived 
from the extragalactic HH region data. 

It has long been known that the D + ^He analysis of Yang et al. (1984) and 
Walker et al. (1991) is likely overly conservative. In both analyses the synthesis 
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Figure 2. The predicted 68% and 95% CL contours for the ^He pri- 
mordial mass fraction with t] constrained by D, ^He and "^Li. Also 
shown is the ztlu range for Yp inferred from the data. 

of new '^He in low mass stars (Iben 1967; Rood 1972; Iben & Truran 1978) 
was neglected. But, Rood, Steigman & Tinsley (1976) had demonstrated that 
such production might dominate the primordial (D + '^He) contribution. Even 
neglecting this contribution, Steigman & Tosi (1992) had followed the evolution 
of ^He in a variety of chemical evolution models and found more ^He survival 
(ffs^ 1/2 rather than 1/4) leading to a higher lower bound to rj. More 
recently, Steigman & Tosi (1995) revisited the "generic" evolution of D and 
^He and, using updated solar system data (Geiss 1993) inferred (for 53 > 1/4) 
^10 ^ 3.1. In a more sophisticated implementation of the "generic" approach, 
Hata et al. (1996) found (for c/3 > 1/4) r/10 > 3.5. 

Although it is still true that rjMiN < Vmax, the increasing lower bound to rj 
increases the lower bound to Ybbn- For rjiQ > 3.1 (3.5), Ybbn > 0.241(0.242). 
An accurate determination of Yp from observations of ^He in low metallicity 
extragalactic HII regions is required for the "^He test. From their analysis of this 
data, Olive & Steigman (1995) derive Yp = 0.232 ± 0.003 where 0.003 is the la 
statistical uncertainty. Thus, at 2a, Y^'^^ < Y^^J^ , faihng the ^He test. It 
should be noted that tjmin and Y^^J^ have already been pushed to their "2(t" 
lower bounds so this discrepancy is at greater than the 95% confidence level. 
The crisis emerges! 

Indeed, Olive & Steigman (1995) used all the data (D, ^He, "^He, ^Li) to 
infer N^, = 2.17 ± 0.27 which deviates from the standard model value (A'^^ = 3) 
by ~ 3a. This crisis for SBBN is reflected in Figure 2 where D, ^He (53 > 1/4) 
and ^Li have been used to bound r], leading to predictions of Ybbn at the 68% 
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Figure 3. The likelihood function (arbitrary normalization) for the 

combined fit (D/'^He, ^He, '''Li) of the data and BBN as a function of 
j^BBN _ At each value of N^^^ the likelihood is maximized for 77. 

and 95% CL. Unless the primordial abundance of ^He has been systematically 
underestimated, the evidence signals a potential crisis for SBBN. 

3. A Statistical Analysis of BBN 

To explore more carefully the consistency of SBBN, my colleagues and I (Hata 
et al. 1995; Thomas et al. 1995) have undertaken the first comprehensive 
statistical analysis of the confrontation between theory and observation. We 
have reexamined the nuclear and weak interactions and their uncertainties and 
have performed a Monte Carlo analysis of the BBN predictions. Indeed, the 
curves in Figure 1 reflect the ilcr uncertainties in the predictions (for Ny = 3) 
of Ybbn, y2P, Utp vs. rj. Prom our Monte Carlos we derive P{A)bbn, the 
probability distributions for the predicted BBN abundances (A). We have also 
reexamined the observational data, accounting for the statistical uncertainties as 
well as attempting to allow for various systematic uncertainties which may arise 
in using the data to infer the distribution, P(A}obs, of primordial abundances. 
These latter uncertainties are not necessarily modelled by gaussian distributions. 
In contrast to previous approaches which treated each element one at a time, we 
may use the information on all light nuclides to form a likelihood function (as a 
function of r] and N^) from P(A)bbn and P{A)obs- The likelihood function, 
maximized with respect to tj at each Ni, is shown in Figure 3 (from Hata et al. 
1995). We derive = 2.1 ± 0.3, consistent with Olive & Steigman (1995). It 
is clear that SBBN {N^y = 3) provides a poor fit to the primordial abundances 
inferred from the data. 
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Figure 4. The 95% CL contours in the vs. r] plane for several 

choices of the ^He survival fraction g3. 

As with the 19th century standard model, SBBN is challenged. As with the 
challenges to the 19th century standard model, there are several options for the 
resolution of this crisis. 

3.1. Is The Chemical Evolution Extrapolation Wrong? 

One source of the challenge to SBBN is the relatively high lower bound to 
7/ imposed by the relatively low primordial abundances of D and ^He inferred 
from solar system and interstellar observations. These stringent upper bounds to 
primordial D and "^He are suggested by many specific chemical evolution models 
(Steigman k. Tosi 1992) as well as the "generic" model for the evolution of D 
and ^He (Steigman & Tosi 1995; Hata et al. 1996). In the latter case, the crisis 
worsens with increasing 53 and/or if stellar production of ^He is allowed for. The 
crisis could be ameliorated if 33 is less than the lower bound (53 > 1/4) adopted 
in the above analyses. In Figure 4 (from Hata et al. 1995), 95% CL contours 
are shown in the vs. rj plane for several choices of 53. If the "effective" 53 
(averaged over stars of all masses and the evolution history of the ISM) is ~ 0.1, 
consistency of SBBN is reestablished. 

3.2. Is The Primordial ^He Abundance Larger? 

An alternate source of the challenge to SBBN is the relatively low abundance of 
primordial ^He inferred from the observations of extragalactic HII regions. By 
allowing only for statistical uncertainties perhaps we've underestimated the true 
uncertainty in Yp. A larger value for Yp could reestablish the consistency of ^He 
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Figure 5. The 95% CL contours in the N^^^ vs. rj plane for several 
choices of the systematic error {AYgys) in the ^He abundance inferred 
from HII region data. 

with D and ^He. Many sources of possible systematic uncertainty in Yp have 
been identified and some have been studied (Davidson Sz Kinman 1985; Pagel 
et al. 1992; Skillman & Kennicutt 1993; Skillman et al. 1994; Copi, Schramm 
Sz Turner 1995; Sasselov & Goldwirth 1995; Olive & Steigman 1995). In Figure 
5 (from Hata et al. 1995) are shown 95% CL contours in the Ni, vs. r) plane for 
several choices of AYgys, where Ybbn = 0.232 it 0.003 + AYgys- If Yp is shifted 
up by ^ 0.010, SBBN may be consistent at the 95% CL. It should, however, be 
emphasized that AYgys may be negative as well as positive; a negative AYgyg 
exacerbates the crisis for SBBN. 

3.3. Is There New Physics? 

By employing N^, as a second parameter, we have allowed for a class of modifica- 
tions of the standard model. If, in addition to three flavors of light, left-handed 
neutrinos {Ni, = 3, SBBN) there are additional light neutrinos or other new 
particles, Ni, > 3 (Steigman, Schramm &; Gunn 1977) and the crisis worsens. 
However, although and arc known to be "light" {rriu < Tbbn)-, accelera- 
tor data on Vj- (ALEPH Collaboration) permits m^^- < 24MeV . As Kawasaki 
et al. (1994) have shown, the presence of a massive, unstable tau neutrino 
with 5 — 10< mi,r < 2AMeV and 0.01< r,yr~ Isec. would correspond to an 
"effective" A^jy < 3. Perhaps the crisis for SBBN is teaching us about extensions 
of the standard model of particle physics. 
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4. Summary and Conclusions 

Primordial nucleosynthesis must have occurred during the early, hot, dense evo- 
lution of a Universe described by the hot big bang model. Therefore, BBN offers 
a test of standard cosmology as well as a probe of particle physics. As with the 
standard model of 19th century physics, over many years SBBN has provided 
support for the standard model of cosmology. Indeed, the success of SBBN in 
predicting the abundances of the light nuclides with only one adjustable param- 
eter rj {Ny = 3) restricted to a narrow range (3< r/io^ 4) while the abundances 
range over some 9 orders of magnitude, is impressive indeed. 

However, as with the 19th century standard model, some clouds have now 
emerged on the horizon. Recent analyses (Copi, Schramm &; Turner 1995; Olive 
&; Steigman 1995; Hata et al. 1995) point to a crisis unless the data are in error, 
or the extrapolations of the data are in error, or there is new physics. Some 
analogies with the crises which confronted the 19th century standard model 
may be instructive. 

(i) Perhaps our extrapolations of the observations of D and ^He from here and 
now to there and then have been naive. Chemical evolution models in which 
more D is cycled through stars and destroyed (without a concommitant over- 
production of ^He) would permit a higher primordial abundance of D, allowing 
a lower rj and Ybbn consistent with Yp. Thus, as with the discovery of Uranus, 
the crisis for SBBN may teach us something new about galactic evolution. 

(ii) Perhaps our estimates of the primordial abundance of ^He are in error be- 
cause we have overlooked some large systematic error in the abundance deter- 
minations. If Yp is larger than the value inferred form the observational data, 77 
may be as large as inferred from D and "'He (see Figure 2) and still Ybbn and 
Yp may be consistent. Then, as with the discovery of Pluto, our crisis may have 
been a false alarm from which, nonetheless, we learn something new. 

(iii) The most exciting possibility, of course, would be that the data are accu- 
rate, the systematic errors small and the extrapolations true. Then, this crisis 
may point us to new physics beyond the standard models of particle physics 
or cosmology. The "window" on a massive, unstable r-neutrino is accessible to 
current accelerators. 

To summarize, then, SBBN [Ny = 3; (73 > 1/4; l^Ygys < 0.005) provides a 
poor fit to the primordial abundances of the light nuclides inferred from current 
observational data (Hata et al. 1995). This crisis is not a cause for alarm, but 
an opportunity to learn something new about astronomy, cosmology, or particle 
physics. 
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